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|~) ' Abstract 

The recurrent nova RS Oph experienced an outburst in 2006, 21 years after its 
previous explosion in 1985, as expected. It was observed at almost all wavelengths, 
^ . and important information about its properties is still being extracted. We present 

lO ! theoretical models of the explosion of this fascinating object, which indicate that the 

JlT^ ' mass of the accreting white dwarf should be very close to the Chandrasekhar mass, to 

^r^ ■ allow for such a short recurrence period. In addition, since models suggest that this 

I — I i nova ejects less mass than it accretes, it is an excellent candidate for a thermonuclear 

^^ ' supernova explosion, in about 10^ — 10^ years from now. We also analyze the emission 



O 



of soft gamma-rays by RS Oph detected with the BAT instrument onboard Swift, 
and with the PCA onboard RXTE. We rule out that this emission has its origin in 
radioactive decays in the expanding nova envelope. 

X' 
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1 Introduction 



Recurrent novae are defined as systems with more than one recorded nova 
outb urst. There are few members of the recurrent novae subclass: just 10 
(see 



Anupama fc Mikolajewskal (Il999l ) and review by iMikolajewskal (120081 ) ) 



divided in two groups according to their orbit al periods. RS Oph belongs to 



the long period group, with Porb = 456 days (jFekel et al.l . l2000l ). which are 
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interpreted as interacting binaries where a hot white dwarf accretes mass from 
a giant secondary; these recurrent novae are a subclass of symbiotic binaries. 
In contrast, classical novae occur in cataclysmic variables, a different type 
of binary system where a white dwarf accretes mass from a main sequence 
companion, overflowing its Roche lobe; typical periods of cataclysmic variables 
are about a few hours. 

RS Oph has undergone various recorded outbursts (1898, 1933, 1958, 1967, 
1985 and 2006, and two additional possible outbursts in 1907 and 1945). In 
the two last outbursts, a lot of observational data was gathered. The recent 
outburst in 2006 offered the opportunity to observe RS Oph at practically all 
wavelengths, radio, IR, optical. X-rays and soft gamma-rays. However, the very 
important UV range was missing, since there was not a satellite "equivalent" to 
lUE flying in 2006. It was a pity, since lU E made very impor tant observations 



in the previous RS Oph eruption in 1985 (jShore et al.l . ll996l ). which could not 
be repeated in 2006. On the other hand, current space missions like Spitzer, 
RXTE, XMM-Newton, Chandra and Swift have provided a wealth of IR and 
X-rays data, of much better quality than what could be obtained at such 
energy ranges in the previous outburst. All in all, there is now an impressive 
data set; its interpretation is a real challenge for theorists and it has provided 
already many exciting results. 

It is worth noting that there is not yet a self consistent model of the RS Oph 
outburst itself, since it is not straightforward to find a set of initial conditions 
leading to a thermonuclear runaway (TNR) with a recurrence period as short 
as 21 years; classical novae recurrence periods are much longer (> 10^ — 10^ 
years), well in agreement with the expected time to reach outburst conditions 
with typical accretion rates and initial white dwarf masses in cataclysmic 
variables. The only known way to reconcile a short interoutburst period -about 
some decades- with a TNR (i.e., explosive hydrogen burning in degenerate 
conditions on top of the white dwarf) , is that a very massive white dwarf, close 
to the Chandrasekhar limit, accretes with a large mass accretion rate. This fact 
gives an additional interest to RS Oph and its relatives, since they are potential 
progenitors of thermonuclear -or type la- supernovae, in a relatively short 
time, if the white dwarf grows in mass in each nova eruption. Our evolutionary 
models indicate that ejected mass is smaller than accreted mass, so confirming 
that possibility. 

In this paper we present our recent attempts to model the basic properties 
of the RS Oph recurrent nova outburst, with special emphasis on the ra- 
tio of ejected to accreted mass, the ejecta chemical composition and the po- 
tential gamma-ray emission. The role played by the initial conditions of the 
white dwarf is underlined. We also discuss the impossibility to explain the 
Swift/BAT hard X-soft gamma-rays detections as related with the contin- 
uum gamma-ray emission coming from electron-positron annihilation in the 



expanding ejecta. 



2 Main observational properties 



The last eruption of RS Oph was recorde d on 2006 February 12.83 UT, reach- 
ing magnitude 4.4 (JNarumi et al.l . 120061 ). Observations at other wavelengths 
have provided important data. For instance, radio observations yielded a new 
detection (not obtained in the previous o utburst) of early emissi on at low radio 
frequencies, clearly of nonthermal origin (JKantharia et al.l . 120071 ) . and high res- 
olution VLB A i mage s, clearly showing an asymmetric expan ding shock wave 
( O'Brien et al.l . l2006l ). Early X-ray observa tions with RXTE ( Sokoloski et al.l . 
20061 ) and Swift/BAT (JBode et al.l . l2006l ) showed hard X-ray emission (up 
to 50 keV), a clear indicator of the blast wave from RS Oph, i.e. an outward 
propagating shock wave consequence of the interaction between the expanding 
ejecta and the red giant wind (the latter much denser than typical circumstel- 
lar material in classical novae). Additional data about the early evolution came 
from IR observations, which also showed the temporal evo lution of the shock 
velocity (JDas. Banerjee fc Ashokl . l2006l : lEvans et al.l . l2007l ). There is a differ- 
ence in the shock velocities deduced from the X-ray and IR data, which has 
been shown to indicate efficient particle acceleration; nonlinear diffusive shock 
acceleration of cosmic rays also explains the decceleration of the blast wave in 
RS Oph, whi ch occured faster than predi cted by the standard adiabatic shock 
wave model (JTatischeff fc Hernang . 120071 ). IR observations also provide an es- 
timate of the ejected mass: ~ 3 x lO^^M© (JDas. Banerjee fc Ashokl . 120061 ). An- 
other important set of observations in X-rays (with Swift /XRT, XMM- Newton 
and Chandra) revealed highly variable soft X-ray emission at about 30 days 
post outburst, and lasting for abou t 60 days, with high temporal variability 



fiBode et all . l2006l : iNess et all . l2007l : iNelson et all . 120081 ). The short duration 



of the su persoft X-ray phase po i nts a. gain to a very massive white dwarf in 
RS Oph (JHachisu. Kato &: Lunai . 120071 ). This supersoft X-ray emission results 
from residual hydrogen burning on top of the white dwarf. 



It is important to know that overabundances of lithium have been found in 
RS Oph in quiescence, and als o in another recurrent no va, T CrB, sharing its 



basic properties with RS Oph fjWallerstein et al.l . 120081 ). 



3 Models 



The model for classical novae explosions is based on the thermonuclear run- 
away ensuing hydrogen degenerate burning, on top of an accreting white dwarf, 
which drives mass ejection at large velocities (hundreds to thousands of km/s). 



For typical initial masses and luminosities of the white dwarf, and mass ac- 
cretion rates from the secondaries, there is an explosion every 10^ — 10^ years. 
Some kind of mixing between the accreted matter and the underlying white 
dwarf should be invoked to explain the large enhancements above solar values 
of CNO nuclei, and also of neon, found in th e ejecta of many n ovae (a list 
of abundances in novae ejecta can be found in lGehrz et alJ ( 19981) ). However , 



no large overabundances of metals were found in RS Oph (jShore et al.l . Il996l ) 



so that not much mixing is expected to occur. The main difficulty with the 
"standard" TNR scenario for RS Oph is the short recurrence period, Prec = 21 
years. Both short Pj-ec and lack of overabundances above solar in the ejecta 
point to fast evolution, only attainable with some particular combinations 
of large white dwarf mass and accretion rate. A first set of nova models on 



white dwarfs with M=1.38 M© were computed by lStarrfield. Sparks fc Truran 



(119851 ) . At ne arly the same epoch, an alternativ e scenario for RS Oph outburst 
was invoked (JLivio. Truran fc Webbinkl . Il986l ). casting doubts on the TNR 
model and assuming a main sequence companion instead of a giant. However, 
it is nowadays accepted that the companion of the white dwarf in RS Oph is 
a red giant and that a TNR is responsible for its outburst. One observational 
fact that proves this scenario is the residual steady H-burning revealed by the 
supersoft X-ray emission detected ~ 30 days after outburst. Other arguments 
come from the shock wave revealed both in X-rays, IR and radio. 



We reanalyse here the TNR model for RS Op h, by means of the hydrody- 
namic code described in iJose fc Hernanzl (119981 ). which follows the evolution 
of the accreting white dwarf from the accretion phase up to mass ejection, 
with complete nucleosynthesis included, and updated nuclear reaction rates. 
The initial conditions are very relevant for the ensuing properties of the ex- 
plosion. An initial mass always larger than 1.2 M© has been adopted, with a 
range of mass accretion rates and white dwarf luminosities. No mixing of any 
kind between the accreted mass (assumed to have solar composition) and the 
white dwarf core has bee n adopted. The main issue i s the i nitial white dwarf 
luminosity; according to lAnupama fc Mikolajewskal (119991 ). white dwarfs in 
symbiotic stars are quite luminuous, as a consequence of steady H-burning 
on their surface. However, such high luminosities do not mean that the whole 
white dwarf core is very hot, but just its outer envelope. Therefore, once the 
purely symbiotic phase ends, a central core temperature-luminosity relation is 
again valid, and relatively low luminosities, such as those required to power a 
TNR, are possible. 



A search for appropriate initial conditions leading to short recurrence periods 
is the first step to model RS Oph. The recurrence period can be written as 
Prcc = AMacc/M = 21 ycars, where AMacc is the accreted mass on top of 
the white dwarf required to power the outburst through a TNR and M is the 
mass accretion rate. The question is which combinations of the initial white 
dwarf mass, M™^, M and the initial luminosity, L^^, lead to a nova outburst 
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Fig. 1. Upper panel: accreted masses required to reach hydrogen ignition condi- 
tions and to power a TNR. Lower paneh recurrence periods for the same set of 
initial conditions as in the upper panel. Asterisks denote the values obtained for 



M = 2 X 10" ^Me/year and L 
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with such a short period of 21 years. We have adopted a range of L™^ between 
10"^ and 1 Lq, and two values for M (2 x 10"'^ and 10"*MQ/year). The initial 
white dwarf mass ranges from 1.25 to 1.38 Mq. The chemical composition of 
the white dwarf should be ONe (and not CO) according to st andard stellar 
evolution (see for instance iRitossa. Garcia-Berro Sz Iberu (1l996l )). 



We show in Figure 1 (upper panel) the accreted masses needed to reach 
hydrogen ignition conditions versus the white dwarf initial mass, for var- 
ious luminosities and M = lO~^M0/year. We also include the results for 
M = 2 X 10~'^MQ/year and L™^ = 1O~^L0. The corresponding recurrence 
periods are shown in the lower panel of the same figure. One sees at first 
glance that the critical accreted mass does not depend only on M™^, as some 
basic approximate formula indicate (since they only require that some pres- 
sure at the base of the accreted envelope -about lo^^"^^ dyn/cm^- is reached 
to power degenerate H-ignition). As expected, an increase in the initial lu- 
minosity reduces the accreted mass needed to power the TNR and thus the 
recurrence period. Another result is that the decrease in accreted mass and 
recurrence period is steeper as the initial white dwarf mass approaches the 
Chandrasekhar limit, since the star becomes easier to destabilize. The main 
issue is that recurrence periods as short as 21 years are really hard to obtain; 
only masses larger than ~ I.SSMq and accretion rates about 10~'^MQ/year 
can provide them. It is worth noting that the initial luminosity should be 
reasonably low, L™^ = 1O~^L0, to guarantee a TNR. Another possibility is a 
smaller accretion rate, lO~^M0/year, with larger L, 1 L0, but the explosion is 
then less powerful and marginally compatible with observations. 

The high accretion rates required agree with some recent predictions from 3D 
simulations, yielding a red giant mass loss rate of about 10~'^Mf7^ /year and an 



accre tion rate of about 10% this value, i.e. 10 ^M0/year (IWalder. Folini fc Shord . 



2008f) . It is worth mentioning, however, that observations by IWorters et al. 
(I2OO7I ) have detected the resumption of optical flickering in RS Oph, indicat- 
ing that mass accretion had been restablished, with deduced accretion rates 
spanning the range between 10~^° and lO~^M0/year, depending on the mech- 
anism of accretion (either wind accretion or Roche lobe overflow); with such 
small accretion rates , recurrence periods \y ould be larger than 21 years by 



orders of magnitude (jJose Sz Hernanzl . Il998l ) . 



In Table 1 we show some properties of the explosion for the most interesting 
cases. The models with M^^=1.35 and 1.38 M©, M = 2 x lO'^Me/year, L^^ = 
1O~^L0 illustrate the effect of the initial white dwarf mass: a slight increase in 
initial mass from 1.35 to 1.38 M© reduces the recurrence period by more than 
a factor of 2 (from 24 to 10 years). The properties of two additional models, 
with M = lO~'^M0/year and two initial luminosities - L™^ = 10~^ and 1 L0 - 
are also displayed, to illustrate the effect of the initial white dwarf luminosity 
and of the mass accretion rate. A decrease in the accretion rate by a factor 
of 20 increases the recurrence period by a similar factor (from ~ 10 to ~ 160 
years), provided that the initial luminosity is unchanged (1O~^L0). On the 
other hand, an increase by a factor of 100 in the initial luminosity (from 10~^ 
to 1 L0) just reduces the recurrence period by a factor of 2.7 (from ~ 160 to 
~ 60 years), provided that the accretion rate and the initial white dwarf mass 
remain unchanged (1.38 M© and lO~^M0/year, respectively). 



Table 1 

Main properties of the RS Oph models: peak temperature, accreted and ejected 
masses, timescales (accretion, recurrence time), total white dwarf mass increase 
and time needed to reach Mchandra- Units are: 10^ K for T, Mq for masses, years 
for times. 
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Peak temperatures (see again Table 1) range around (2 — 3) x 10^ K, which 
are moderate but similar to those of classical novae; in contrast, the total 
time elapsed since accretion starts up to ignition is much shorter, as already 
emphasized. The corresponding accreted and ejected masses are quite small, 
~ IQ-® Mq or even smaller (~ 10^^ Mq for the model with large initial 
luminosity, 1 Lq). However, the case with large initial luminosity does not 
give a very energetic explosion and would not well represent RS Oph (but it 
is not too bad concerning the recurrence period, 58 years, just a factor of ~ 3 
above the observed value). The two models with M = 2 x lO^^MQ/year are 
thus those that best reproduce RS Oph main observed features. The maximum 
bolometric luminosities reach ~ (3 — 5) x IO^Lq, the average global kinetic 
energy is about 10"^^ erg, with average velocities of (3 — 4) x 10^ km/s, and 
ejected masses are ~ 10"^ Mq, which indicate that this nova is a bit less 
energetic than classical ones, and with a smaller content of radioactive nuclei 
(as shown below). 

An important result is that the accreted masses are larger than the ejected 
ones, so that the mass of the white dwarf increases: Maec- Mejee ranges between 
0.7 and 1.7 x 10^® Mq, leading to an expected time to reach the Chandrasekhar 
mass (Mchandra= 1-4 Mq) ranging between 3 and 7 xlO^ years for the two 
models with M = 2 x IQ-^MQ/year. The other models would need much 
larger times to reach the Chandrasekhar mass, and thus would not be viable 
type la progenitors. 

The nucleosynthesis in our models provides large hydrogen and helium abun- 
dances in the ejecta (mass fractions ~ 0.6 and ~ 0.4, respectively). There is 
not large overproduction with respect to solar of metals, in agreement with 
observations, and contrary to what is obtained in general for classical no- 



vae. The main reason for such a difference is the lack of initial mixing and 
the rapid evolution to the runaway. Regarding ^Li, it is not at all overpro- 
duced. The key issue in the synthesis of ^Li is the amount of ^Be (which is 
transform ed into ^Li in t h e ejec t a after an electron capt ure) that survives the 



TNR (see lHernanz et al.l (Il996l ): iJose fc Hernanzl (119981 )). At the early stages 



of the TNR, ^Be is efficiently destroyed by proton capture reactions, but as 
soon as enough T is reached (about (5 — 7) x lO^K) photodisintegration of 
^B drives a pseudo-equilibrium between ^Be and ^B. Hence, the final amount 
of ^Li strongly depends on the level of ^Be destruction prior to this pseudo- 
equilibrium. The faster the TNR, the shorter to reach the pseudo-equilibrium 
stage, and hence, the larger the final amount of ^Li. This explains why CO 
models lead in general to larger amounts of ^Li in the ejecta as compared with 
ONe models (since mixing with the underlying CO WD brings ^^C fuel into 
the envelope which speeds up the TNR dramatically). Also, the increase of 
strength in explosions hosting massive white dwarfs explains why very mas- 
sive ONe white dwarfs (close to the Chandrasekhar mass) eject larger amounts 
of ^Li than not so massive ONe white dwarfs, but never reaching solar val- 
ues or larger. Thus, the Li overabundance found in RS Oph in quiescence 
( IWallerstein et al.l . l2008l ) is hkely not synthesized during the nova outburst. 



There are small amounts of t he radioactive nuclei relevant for 7-ray emission 
(^^N, ^^F, ^Be and ^^Na, see iGomez-Gomar et al.l ( 19981 )). i.e., much smaller 
yields in mass fraction (and also absolute yields, because of the small ejected 
masses) than for classical novae. In addition, the envelope becomes trans- 
parent to 7-rays later than in normal novae, because of the smaller global 
kinetic energy. All in all, the expected gamma-ray fluxes are very small. This 
rules out the interpretation of the observed soft 7-rays from RS O ph, with 



the BAT instrument on board the Swift satellite (JBode et al.l . 120061 ). as due 



to the decay of radioactive nuclei. Shocks within the e jecta should provide a 



2006 



better explanatio n of the observed hard X-ray emission (iSokoloski et al.l 
Bode et al.l . l2006l ). In fact, one could know beforehand that Swift /BAT obser- 
vations of RS Oph can not be interpreted as due to ra dioactive nuclei decay, 



by two reasons: the spectral shape and the light curve (jGomez-Gomar et al. 



19981 : iHernanz. Gomez-Gomar fc Josel . |2002| ) . The shapes of the 7-ray spectra 



expected from radioactive decay at various epochs are such that there is less 
emission in the 15-25 keV energy band than in the 25-50 keV band, due to 
the sharp spectral cutoff at low energies produced by photoelectric absorption. 
The co ntrary was observe d with Swift/BAT and RXTE, as shown in Figure 
1 from iBode et al.l (120061 ). Concerning the 7-ray light curve, our models in- 
dicate that the maximum in visual flux occurs later than the maximum in 
7-rays (simultaneous to peak temperature), for all types of novae, again in 
contradiction with the observations with Swift/BAT, where the peak in hard 
X/soft 7-ray emis sion occurred after the optical outburst (see Figure 1 from 



Bode et all (l2006f )) 



4 Discussion and conclusions 



Models of the RS Opli explosion should account for its main observed proper- 
ties, e.g., very short recurrence period (21 years), small ejected mass, roughly 
solar chemical composition of the ejecta and maximum luminosities and ejecta 
velocities similar to those of classical novae. Hydrodynamic models of massive 
(larger than ~ 1.35 M©) accreting white dwarfs with large accretion rates - 
ranging from M = 10~* to 2 x lO~'^M0/year - well reproduce these features. 
But there are still some aspects to be better understood, as the need of ini- 
tially low white dwarf luminosities, as compared with some observations at 
quiescence. We want to point out a few important puzzles. First of all, very 
massive white dwarfs should be made of ONe (if they were born as massive 
white dwarfs); therefore, they pose a problem as scenarios for type la su- 
pernovae, because ONe white dwarf s rather collapse to a neutron star than 



explode, when reaching Mchandra (see iGutierrez. Canal &: Garcia-Berrd (120051 ) 



and references therein). One could envision that the white dwarf was initially 
born as a CO white dwarf, but then it should have been less massive than 
1.1 M0 and be able to reach 1.35-1.38 M© through successive recurrent nova 
eruptions. However, this seems unrealistic because models do not reproduce 
such behaviour. Second, the mass left after the recurrent nova explosion, i.e.. 
Mace- Mejec (scc abovc) should be compatible with the duration of the sup er 



soft X-ray emission phase (around 60 days, seelHachisu. Kato fc Lunal (120071 )) 



According to our models (jSala fc Hernana . l2005l ) . the remnant envelope mass 



for a hydrogen mass fraction of 0.6 and mass of the white dwarf larger than 1.3 
Mq should be around 10"''' Mq, and surely smaller than 10~^ Mq. Therefore, 
the mass of the white dwarf is more plausibly 1.38 than 1.35 M© (see Table 
1). Finally, the hard X/soft 7-rays detected with Swift /BAT and RXTE in 
the early phases of the explosion could not come from electron-positron an- 
nihilation emission (with positrons coming from the decay of the /^"'"-unstable 
nuclei ^^N, ^^F in the ejecta), both because of the tiny amount of radioactive 
nuclei, the too late appearence of such emission and its spectral shape. Deeper 
studies of both the thermonuclear runaway phase and the interaction between 
the ejecta and the dense stellar wind of the red giant are in progress, to better 
understand the wealth of observational data available for the fascinating RS 
Oph recurrent nova. 



Acknowledgements 



This work has been supported by the Spanish MEC and by the Catalan DURSI 
through grants ESP2007-61593, AYA2007-66256 and 2005SGR00378, and by 
FEDER funds. We thank Vincent Tatischeff for many valuable discussions. 



References 

Anupama, G.C. & Mikolajewska, J. 1999, A&A, 344, 177 

Bode, M.F. et al. 2006, ApJ, 652, 629 

Das, R., Banerjee, DR.K. & Ashok, N.M. 2006, ApJ, 653, L141 

Evans, A. et al. 2007, MNRAS, 374, LI 

Fekel, F.C., Joyce, R.R., Hinkle, K.H. & Skrutskie, M.F. 2000, Astron. J., 119, 

1375 
Gehrz, R.D., Truran, J.W., Williams, R.E., Starrfield, S. 1998, PASP, 110, 3 
Gomez-Gomar, J., Hernanz, M., Jose, J. & Isern, J. 1998, MNRAS, 296, 913 
Gutierrez, J., Canal, R. & Garcia-Berro, E. 2005, A&A, 435, 231 
Hachisu, I., Kato, M. & Luna, G.J.M. 2007, ApJ, 659, L153 
Hernanz, M., Jose, J., Coc, A. & Isern, J. 1996, ApJ, 465, L27 
Hernanz, M., Gomez-Gomar, J. & Jose, J. 2002, New Astr. Rev., 46, 559 
Jose, J. & Hernanz, M. 1998, ApJ, 494, 680 
Kantharia, N.G., Anupama, G.C., Prabhu, T.P., Ramya, S., Bode, M.F., 

Eyres, S.P.S. & O'Brien, T.J. 2007, ApJ 667, L171 
Livio, M., Truran, J.W. & Webbink, R.F. 1986, ApJ, 308, 736 
Mikolajewska, J. 2008, in Proceedings of the meeting: RS Ophiuchi (2006), ed. 

by N. Evans, M. Bode and T. O'Brien, Astron. Soc. of Pacific Conference 

Series, in press (astro-ph/0803.3685) 
Narumi, H. et al. 2006, lAU Circular 8671 
Nelson, T., Orio, M., Cassinelli, J. P., Still, M., Leibowitz, E. & Mucciarelli, 

P. 2008, ApJ, 673, 1067 
Ness, J.-U. et al. 2007, ApJ, 665, 1334 
O'Brien, T.J. et al. 2006, Nature, 442, 279 
Ritossa, C, Garcia-Berro, E. & Iben, L 1996, ApJ, 460, 489 
Sala, G. & Hernanz, M. 2005, A&A, 439, 1061 

Shore, S.N., Kenyon, S.J., Starrfield, S. & Sonneborn, G. 1996, ApJ, 456, 717 
Sokoloski, J. L., Luna, G.J.M., Mukai, K. & Kenyon, S.J. 2006, Nature, 442, 

276 
Starrfield, S., Sparks, W.M. & Truran, J.W. 1985, ApJ, 291, 136 
Tatischeff, V. & Hernanz, M. 2007, ApJ, 663, LlOl 

Wallerstein, G., Harrison, T., Munari, U., Vanture, A. 2008, PASP 120, 492 
Walder, R., Folini, D. & Shore, S.N. 2008, A&A, 484, L9 
Worters, H.L., Eyres, S.P.S, Bromage, G.E. & Osborne, J.P. 2007, MNRAS, 

379, 1557 



10 



